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I. ABSTRUCT 


The inconsistency about the degree of geometrical frustration has been a long issue in 
AV2O4 (A = Zn, Cd and Mg) compounds, which arises from the two experimental results: 
(i) frustration indices and (ii) magnetic moments. In the present study, we try to understand 
such inconsistency by using ab initio electronic structure calculations. The orbital degrees of 
freedom are found to play an important role in understanding the geometrically frustrated 
magnetic behaviour of these compounds. The magnitude of the maximum calculated values 
of orbital magnetic moment per formula unit for ZnV204, MgV204 and CdV204 compounds 
are found to be ~ 1.54 /i^, ~ 0.92 /i^ and ~ 1.74 /i^, respectively. The inclusion of the orbital 
and spin angular momenta for calculating the frustration indices improves the understanding 
about the degree of geometrical frustration in these compounds. The calculated values of the 
frustration indices (/j) are largest for MgV204 and smallest for CdV204 for 3 . 3 < U < 5.3 eV. 
In this range of U, the calculated values of AM2=Mtotai-Mexp (where, Mtotai=Mspin-|Morbitai|) 
are also found to be largest for MgV204 and smallest for CdV204. Hence, the consistency 
about the degree of geometrical frustration, which arises from the / j as well as from the AM2 
is achieved and improves the understanding about the degree of geometrical frustration in 
these compounds. Calculated values of band gap in this range of U are found to be closer to 
that of experimentally observed values for all three compounds. The absolute values of the 
nearest neighbour exchange coupling constant (Jnn) between V spins are found to be largest 
for MgV204 and smallest for CdV204, which indicate that the calculated absolute values of 
the Curie-Weiss temperature { 0 cw)j are highest for MgV204 and smallest for CdV204 for 
3 . 3 < U < 5.3 eV. In this range of U, the magnetic transition temperature (Tn)j is found 
to be ~150 K, ~60 K and ~22 K for MgV204, ZnV204 and CdV204, respectively, which 
shows that the order of (Ttv) j is similar to that of (TAr)exp for these compounds. Hence, all 
the magnetic properties studied in the present work are well explained in these spinels for 
3 . 3 < U < 5.3 eV. This work is expected to provide a valuable input in understanding the 
geometrically frustrated magnetic behaviour for those systems for which the orbital part of 
the angular momenta are not quenched. 
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II. INTRODUCTION 


The ground state properties in the condensed matter physics have been well described 
from long time by the density functional theory (DFT) approach.^®^ At least in the metallic 
systems, such properties are described well either by local density approximation (LDA) or 
generalized gradient approximation (GGA) based on the DFT. However, it is well known that 
these functionals underestimate the orbital moments for strongly correlated systems, which 
are induced by the strong spin orbit coupling (SOG).®^This can be resolved by adding or¬ 
bital dependent Hartree-Fock (HF) potential to the LDA/GGA, so called LDA-I-f//LDA-|-[/ 
approximation.f^HIo] \ major problem in this approach is that the electron-electron interac¬ 
tion has already been included in LDA/GGA potential. Hence, the simple addition of HF 
potential leads to double counting (DG). Here, the best way is to identify the mean-field-part 
of the HF potential and subtract it, leaving the orbital dependent correction to the mean 
held type LDA potential. Gzyzyk and Sawatzky suggested a scheme that is true for uniform 
occupancies, so called around mean held (AMF) and is applicable to the weakly correlated 
systems.^ For strongly correlated systems, AMF DG scheme is not valid due to the lack of 
uniform occupancies. For these systems, one can prefer the fully-localized (FL) DG scheme, 
where the average ehect for a localized state is subtracted with integer occupation number.® 

Geometrical frustration is always a challenging problem in the strongly correlated sys¬ 
tems. In these systems it arises due to the interactions between spin degrees of freedom in 
a lattice, which are incompatible to that of the essential crystal geometry. In the highly 
geometrically frustrated magnets, frustration suppresses the long-range magnetic order and 
leads to a degenerate manifold of ground states. Degeneracy in the frustrated magnetic sys¬ 
tems yield diherent complex ordering structures, spin liquid states and spin ice states.^^nSl 
Gharge ordering phenomena are also much ahected by the geometrical frustration.^ Af¬ 
ter 1980 ’s, spin systems on the pyrochlore lattice (an example of geometrically frustrated 
structures) have been studied in more detaih^^^UMlQ] classical Heisenberg spin sys¬ 

tems, a very strong geometrical frustration is anticipated due to the antiferromagnetically 
coupled spins on a pyrochlore lattice, which do not show long-range order at any nonzero 
temperature 

Vanadium spinels, AV2O4 (A = Zn, Gd and Mg) with the face-centered-cubic structure 
at room temperature are an interesting example of geometrically frustrated systems, which 
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show a variety of physical properties.^^^^ In these systems a pyrochlore lattice is formed due 
to the corner sharing network of tetrahedra with magnetically coupled V atoms that gives 
the geometrical frustration.®^ As opposed to the other geometrically frustrated systems, 
vanadium spinels show the long range antiferromagnetic ordering at low temperature.® In 
these compounds, orbital ordering is found to be responsible for the structural transition 
that leads to a long range magnetic ordering due to the lifting of geometrical frustration 
in the tetragonal structure.®® The structural transition temperature (T5) for CdV204 
(~97 K)>MgV204 (~65 K)>ZnV204 (~50 K) compound.®®®®® However, the mag¬ 
netic transition temperature (T^r) for MgV204 (~42 K)>ZnV204 (~40 K)>CdV204 (~35 
K) compound.®®®®® The values of Tat for these compounds are always found to be less 
than Ts- However, it is expected that the values of Ts and should be same for these 
compounds because the geometrical frustration is removed by the above proposed mecha¬ 
nism. The fact that Tp^<Ts in these compounds, indicates the presence of certain degree 
of geometrical frustration in the tetragonal phase of the compound. 


The presence of geometrical frustration in the above mentioned compounds can be at¬ 
tributed to the following two experimental facts. First one is the frustration index, dehned 
s-s f s=\G>cw\/TN, where |6>cm| and Tat are the Curie-Weiss temperature and the magnetic 
ordering temperature, respectively. Experimentally reported values of Ocw for ZnV204, 
MgV204 and CdV204 are ~850 K ~600 K and ~400 K, respectively.®® However, in some 
literature these values are reported to be different for these compounds depending upon 
the experimental conditions.®®® In general, one can use the parameter as a measure 
of strength of geometrical frustration in the system.®®® The values of fs for ZnV204, 
CdV204 and MgV204 are found to be 21 . 3 , 11.4 and 14 . 3 , respectively.® Among these com¬ 
pounds, ZnV204 is largest frustrated and CdV204 is least frustrated as large value of fs 
indicates the more frustration in the system. The second one is based on the difference be¬ 
tween the spin magnetic moment per formula unit (MM/f.u.) [Mspin] and the experimentally 
observed MM/f.u. (Mgxp), which is denoted by AMi=Mspin-Mexp in the antiferromagnetic 
phase of these compounds. Mgpin for these compounds is found to be 4 However, Mgxp 
for CdV204, ZnV204 and MgV204 compounds are reported to be 2.38 /is, 1.26 and 
0.94 hb-i respectively.®®®® Large values of AMi is explained in two different ways. One 
way of understanding is to attribute this difference to the activeness of geometrical frus¬ 
tration in these compounds, which results in the quantum fluctuations responsible for the 
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reduction of Mexp drastically as reported in MgV204.l^ If this is the case, then the degree 
of the geometrical frustration for MgV204>ZnV204>CdV204. Also, the order of geometri¬ 
cal frustration that arises from the experimentally reported frustration index is found to be 
ZnV204>MgV204>CdV204 as described in the reference [ 43 ]. Hence, both the results show 
the inconsistency about the degree of geometrical frustration in these compounds. Another 
way of understanding will be by considering contribution of the magnitude of orbital part 
of MM/f.u. (denoted by |Morbitai|) reported by Maitra et al. for ZnV204.l^ The calculated 
values of the |Morbitai| for ZnV204, MgV204 and CdV204 compounds are found to be ~ 1.50 
fiB, ~ 0.4 fiB and ~ 0.4 respectivelyl ^^ * ^^ * ^^ Even by considering the iMorbitail, the total 
calculated MM/f.u. (Mtotai) is still larger than the Mexp, indicating that the certain degree 
of geometrical frustration is still present in these compounds. The inclusion of |Morbitai| does 
not improve the situation. 


The above discussion clearly suggests that the current understanding of these compounds 
is not sufficient to resolve the contradiction about the degree of geometrical frustration 
present especially in ZnV204 and MgV204. In the present work we have tried to resolve this 
contradiction by carrying out detailed LDA-|-[/-|-SOC calculations. Our results clearly show 
the importance of orbital degrees of freedom along with the spin in understanding the mag¬ 
netic properties of these compounds. The inclusion of orbital and spin angular momenta for 
calculating the frustration indices improves the understanding about degree of geometrical 
frustration present in these compounds. The calculated values of the frustration indices fj 
thus obtained for 3 . 3 < U < 5.3 eV provide the largest value for MgV204 and smallest for 
CdV204. This behaviour is similar to that expected from the AM2. In this range of U, 
the magnitude of Curie-Weiss temperature [{Ocw)j] are highest for MgV204 and smallest 
for CdV204, which indicate that the magnitude of the nearest neighbour exchange coupling 
constant (Jnn) among V spins are largest for MgV204 and smallest for CdV204 compounds. 
Also, the order of magnetic transition temperature [(Ttv) j] is found to be consistent with the 
experimentally reported order of magnetic transition temperature (TTv)exp for these spinels 
for above mentioned range of U. 
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III. COMPUTATIONAL DETAILS 


In present work, the ferromagnetic (FM) and antiferromagnetic (AFM) electronic- 
structure calculations of AV2O4 (A = Zn, Cd and Mg) compounds are carried out by using 
the state-of-the-art full-potential linearized augmented plane wave (FP-LAPW) method.^ 
The atomic positions and lattice parameters used in the calculations for every compounds 
are taken from the literature.®*^^® All these calculations are performed in the tetragonal 
phase for which Perdew -Wang/Ceperley -Alder exchange correlation functional has been 
used.l^ In order to calculate the nearest neighbour exchange coupling constant {Jnn), we 
have considered the two magnetic orderings: FM and AFM (not the experimentally ob¬ 
served structure) ordering of the spins of the four V atoms in a primitive unit cell. For 
AFM ordering, the primitive unit cell consists of two up and two down spins on the four V 
atoms. The effect of on-site Coulomb interaction among V 3 d electrons is considered within 
LDA-|-f/ formulation of the density functional theory.^ Normally in this method U and J 
are used as parameters. However in our calculations, only U is used as a free parameter and 
the value of J is calculated self-consistently as described in the reference [ 7 ]. FL and AMF 
DC schemes have been used in these calculations.® Both the DC schemes give the similar 
results. Thus in order to avoid the repetition of the results, we have discussed the results 
corresponding to the FL DC scheme in the rest of the manuscript. The muffin-tin sphere 
radii used in the present work are 2 . 0 , 2 . 46 , 1 . 39 , 2.0 and 1.54 Bohr for Zn, Cd, Mg, V and 
O, respectively. (6,6,6) k-point mesh size is used in the calculations. In order study the role 
of orbital degrees of freedom, SOC is included self-consistently in LSDA-|-f/ calculations, 
where no external parameters have been used. The effect of the spin canting has been found 
to be small for these spinels both theoretically and experimentally.^^^®^!] Hence, collinear 
calculations have been performed here by varying U from 2-6 eV, where the direction of the 
magnetization is set along the z-axis. In this range of U, every calculations are started from 
the different combinations for converged electron densities and potentials corresponding to 
different values of U for every compounds. The convergence target of total energy has been 
achieved below 10 “^ Hartrees/cell. 
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IV. RESULTS AND DISCUSSIONS 


The crystal structure and splitting of d levels of V ion of AV2O4 (A = Zn, Cd and Mg) 
compounds in presence of octahedral crystal held and SOC are shown in the Fig. l(a and 
b). It is clear from the Fig. 1 (a) that each V ion is located at the center of the edge sharing 
octahedra which create a crystal held. Such a crystal held splits the d levels into lower 
energy t2g levels and higher energy eg levels in the cubic phase as shown in the Fig. 1 (b). 
Because of ignoring the small trigonal distortion present in these systems, these t 2 g levels 
are normally assumed to be threefold degenerate. The t2g levels are generally represented by 
orbital angular momentum 1=1. The degeneracy of t 2 g levels is further lifted by the SOC, 
where the lower value of total angular momentum j=l/2 corresponds to the lower energy 
state. 

In order to see the ehect of SOC for all the three compounds, we have hrst calculated the 
SOC energy (the energy diherence between LSDA+U+SOC and LSDA+U energies for these 
compounds) roughly. Here, as a representative, we have shown the calculated value of SOC 
energy only for U =4 eV. The magnitude of rough estimated SOC energy per formula unit 
for ZnV204, CdV204 and MgV204 are ~ 0.5 eV, ~ 0.5 eV and ~ 0.2 eV, respectively for U =4 
eV. Now, in order to compare the calculated MM to the experimentally observed MM, we 
have used the FM structure instead of experimentally observed AFM structure in the present 
calculations. This is because of the following reasons. It is important to note that the various 
physical quantities {fs, Ocw and T^r) for these compounds are observed experimentally for 
spin, S'=l. Now in order to compare our results with the experimental data (discussed later 
of the manuscript), FM structure calculations are best for calculating the MM as compared 
to the experimentally observed AFM structure calculations for these compounds. This is 
because of the fact that in AFM calculations, the total MM in the interstitial region is 
almost zero and hence the Mgpm (comes only inside the muffin-tin spheres) is always less 
than 4.0 hb for U=2-Q eV, which indicates that S < 1 . However, in FM calculations values 
of Mgpin comes out to be 4.0 fis (both inside the muffin-tin spheres as well as from the 
interstitial region) as per expectation for whole range of U, which indicates that the V atom 
is in 3 -|- ionic state with spin, S=l. However, the orbital MM is calculated only inside the 
muffin-tin spheres because there is no way of calculating the orbital MM in the interstitial 
region. It is important to note that with increase in the muffin-tin sphere radius of V atom 
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more than that used in the present calculations, only a slight changes in the orbital MM 
has been observed, which indicates that the contribution of orbital part of MM is almost 
negligible from the interstitial region. Also for calculating Mtotai (comes from both spin and 
orbital part of MM), we have also calculated the |Morbitai|- 

The plot of MM/f.u. as a function of U is shown in the Fig. 2 (a-c). It is evident 
from the hgure that the values of the Morbitai are negative for these compounds, which 
indicate that the direction of the orbital MM is opposite to the spin MM consistent with 
the Hand’s rule for less than half-hlled cases. The behavior of MM/f.u. for ZnV204 is 
similar (different) to that of MgV204 (CdV204) for whole range of U studied here. The 
|Morbitai| shows a small dependence on U from 2 to 3 eV for both ZnV204 and MgV204 
compounds. In this range of U both compounds gives almost equal values of | Morbitai I, 
where it increases from ~ 0.10 to ~ 0.20 ^b- However, it increases sharply to ~ 1.52 ( 0 . 78 ) 

for ZnV204 (MgV204) compound at 77 = 3.5 eV. The value of |Morbitai| remains almost 
same for ZnV204 and increases to ~ 0.92 hb for MgV204 as U changes from 3.5 to 4 eV. For 
CdV204, it increases from ~ 0.82 to ~ 1.74 hb as U varies from 2 to 3.5 eV. For U >4 eV, it 
decreases continuously upto ~ 0.68 /i^, ~ 0.42 hb and ~ 0.62 /i^ for ZnV204, MgV204 and 
CdV204, respectively at f/=6 eV. The possible cause of changes of |Morbitai| by varying U 
observed in these compounds is discussed below as: In regular octahedron symmetry hve fold 
degenerate d orbitals split into lower energy triply degenerate t2g levels and higher energy 
doubly degenerate e^ levels. For large crystal held splitting, the lower energy t2g levels are 
normally represented by angular momentum 1=1 (/’^=-l, 0, +1) and higher energy e^ levels 
by (/^=- 4 , + 4 ). In tetragonal distortion with c/a < 1 , t2g levels are split into lower 
energy dxy (F=0) orbital and higher energy degenerate dxz and dyz (/^=-l, +1) orbitals. 
Similarly, e^ levels split into dx^-y^ and d^^ orbitals.®^^® Thus the total orbital angular 
momentum (L) of V ions are decided by the number of electrons in each of these hve d 
orbitals. For example, if we consider pure ionic picture along with full orbital polarization 
for these compounds, then one out of two V 3 d electrons occupies the dxy orbital and second 
electron occupies either dxz or dyz orbital depending upon the site as these compounds show 
orbital ordering. Hence, result L =1 for d electrons. However, presence of certain degree of 
hybridization breaks the fully ionic picture of the compounds and 3 d electrons may occupy 
all the hve d orbitals depending upon crystal held splitting energy, Hunds coupling energy, 
on-site Coulomb interaction, etc. These electrons occupancies of various orbitals decide the 




L of the V ions. This discussion clearly suggests that the localization of electrons (which 
changes the total number of d electrons at V sites) along with their distribution in various 
d orbitals finally decide the values of L for V ion. Hence, with increasing U from 2 to 4 
eV, the total number of d electrons (means fractional occupation of all five d orbitals) of 
V ion increases, which result the increase in |Morbitai| because of the more localization of 
the electrons as per expectation.®^^ Above U =4 eV, one can expect that the decrease in 
|Morbitai| is due to the delocalization of electrons. But this is not the case. This is because for 
U >4 eV, the total number of d electrons of V ion remains constant upto U =6 eV. However, 
the redistribution of electrons in all hve d orbitals above U=A eV have been observed in 
the density matrix, which is expected to be responsible for decrease in |Morbitai|- This type 
of non-monotonic behavior of |Morbitai| with increasing U from 2 to 6 eV is expected in 
these compounds as shown by Bultmark et al. in ferromagnetic PuP and US systems.^ The 
maximum calculated value of |Morbitai| for ZnV204 is found to be ~ 1.54 /i^ at 17=4 eV, 
which is in good agreement with that calculated by Maitra et al. as shown in the Table I.®' 
However, it is found to be ~ 0.92 at 17=4 eV and ~ 1.74 /i^ at 17 = 3.5 eV for MgV204 and 
CdV204, respectively and is larger than that predicted by different groups as also shown in 
the Table i.EzEsISs] value of |Morbitai| depends on the exchange correctional functional, 
on-site exchange parameter J, muffin-tin sphere radii and type of calculations (collinear 
or non-collinear), which are different in the present calculations as compared to the earlier 
calculations. From above discussion, it is clear that the calculated values of |Morbitai| are 
different for these compounds for whole range of U studied here. This is due to the fact that 
the distorted VOg octahedron contains V -0 bonds and V-O-V angles, which are different 
for these compounds and are partly attributed to the difference in the ionic radii at the 
A=Cd^’'', Mg^+ and Zn^+ site. Therefore, the values of |Morbitai| are expected to be different 
for all three compounds because it comes from the distribution of electrons among the V 
3 cl orbitals. Such a distribution of electrons in these orbitals are greatly influenced by the 
V -0 bonds and V-O-V angles. It is evident from the hgure that the contribution of the 
|Morbitai| to the Mtotai (where, Mtotai=Mspin-|Morbitai|) is more for CdV204 as compared to 
ZnV204 and MgV204 compounds for 17 = 2-4 eV. Above f /=4 eV, its contribution to the 
Mtotai is largest for ZnV204 and smallest for MgV204 compound. In these compounds, the 
values of M^tai show the similar behavior to that of Morbitai for whole range of U. For U <3 
eV, it is almost equal to the Mgpin for ZnV204 and MgV204 compounds, which is due to 
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the small contribution of |Morbitai|- However, for CdV204, it decreases from ~ 3.18 to ~ 2.26 
as U changes from 2 - 3.5 eV. For ZnV204 (MgV204) compound, it decreases sharply to 
~ 2.48 (~ 3 . 22 ) /is at 17 = 3.5 eV. The value of |Mtotai| remains almost same for ZnV204 and 
decreases to ~ 3.08 ps for MgV204 as U changes from 3.5 to 4 eV. For U >4 eV, it increases 
continuously upto ~ 3.32 fj,B, ~ 3.58 ps and ~ 3.38 fiB for ZnV204, MgV204 and CdV204, 
respectively at U =6 eV. Here, it is interesting to compare the calculated values of Mtotai to 
that of experimentally observed values of Mexp for these compounds, which are also shown 
in the Table 1 . The value of Mgxp has been taken from the references [ 30 , 33 , 34 , 41 ]. The 
calculated values of Mtotai is larger than Mgxp for ZnV204 ( 1.26 and MgV204 ( 0.94 /i^) 
compounds for all values of U. However, for CdV204 compound, the calculated values of 
Mtotai is larger than the Mgxp ( 2.38 hb) for U <3 eV and U >4 eV. For 3 < U <4 eV, the 
values of Mtotai is almost equal to Mgxp. The value of Mexp<Mtotai suggests the presence of 
certain degree of geometrical frustration in these compounds. 

In order to make this point more clear, we have plotted the difference between the Mtotai 
and Mgxp (AM2=Mtotai-Mgxp) in the Fig. 3 . It is clear from the hgure that the behavior 
of AM2 for both ZnV204 and MgV204 compounds are almost similar for whole range of 
U. For ZnV204 (MgV204), AM2 decreases slowly from ~ 2.64 (~ 2 . 96 ) to ~ 2.54 (~ 2 . 88 ) 
Pb when U changes from 2 to 3 eV. At U= 3.5 eV, it decreases sharply to ~ 1.22 (~ 2 . 28 ) 
/is for ZnV204 (MgV204) and remains almost constant upto 17=4 eV. Above 17=4 eV, 
it increases continuously upto ~ 2.06 (~ 2 . 64 ) fiB for ZnV204 (MgV204) at 17=6 eV. For 
CdV204, AM2 decreases continuously from ~ 0.80 to ~ 0.0 fiB as U varies from 2 to 4 eV 
and then increases upto ~ 1.0 fiB at 17=6 eV. It is also clear from the figure that the AM2 
for MgV204>ZnV204>CdV204 for all values of U. Here, we attribute AM2 to the degree 
of activeness of geometrical frustration in these compounds which results in the quantum 
fluctuations responsible for the reduction of MM. Hence among these compounds, the degree 
of geometrical frustration is largest for MgV204 and smallest for CdV204. From above 
discussion, it is clear that the order of AM2 is similar as that of AMi for these compounds, 
which indicates that even by including the Mgrbitai to the Mtotai does not change the order of 
degree of geometrical frustration. Hence, the degree of geometrical frustration, which arises 
from AM2 is inconsistent to that observed from the experimentally reported fs- 

It is important to note that the experimentally reported values of fs in these compounds 
are calculated by the ratio between the Curie-Weiss temperature and the antiferromagnetic 
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transition temperatnre as, 


fs — \{^Cw)exp\/{TN)exp (1) 

The valnes of { 0 cw)exp are calcnlated (withont inclnding the orbital angnlar momentnm 
L) by the following formnla,^ 


(OcwU, = + 1 )A (2) 

ok^v 

Where N/V, gs, S and A are the nnmber of magnetic atoms per nnit volume, 

Lande g factor for the total spin angular momentum, Bohr magneton, Boltzmann constant, 
total spin angular momentum and Weiss molecular held constant, respectively. Substituting 
Eqn. ( 2 ) in Eqn. ( 1 ), we get 


fs 


8iV/i| , 

SksV {T]S[)exp 


( 3 ) 


for 5=1 and gs=‘^- 


For all three compounds, we have seen above that the contribution from the orbital part 
of MM is signihcant and can not be neglected. So, the above approximation is not expected 
to give correct values of the frustration indices. The best way for calculating the frustration 
indices is to include the contribution from both orbital and spin angular momenta. Now in 
order to calculate the frustration indices for these compounds, we have replaced the total 
spin angular momentum S by the total angular momentum J in Eqn. ( 2 ) and after using 
Eqn. ( 1 ), we get 


fj 


N{gjg,B)‘^ 

3 k bV (T/v)exp 


J(J+ 1 )A 


( 4 ) 


Where fj, J=\L-S\, L and gj=l+ frustration index (by in¬ 

cluding both orbital and spin angular momenta), total angular momentum for less than half 
hlled d orbitals, total orbital angular momentum and Lande g factor for the total angular 
momentum, respectively. Dividing Eqn. ( 4 ) by Eqn. ( 3 ), we get 


fj = 


gjJ{J + l) 


fs 


( 5 ) 


The experimentally observed values of fs are given in the Table 1 . 

Now, using Eqn. ( 5 ), we have calculated the values of /j by including both orbital and 
spin angular momenta for these compounds by varying U from 2 to 6 eV. In Fig. 4 , we have 
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plotted the f j for these compounds as a function of U. It is evident from the hgure that the 
values of /j are largest for ZnV204 and smallest for CdV204, for U < 3.3 eV. However, for 
3.3 <U < 5.3 eV, the order of fj is MgV204>ZnV204>CdV204. For U > 5.3 eV, the order 
of /J in these compounds is similar to that observed for U < 3.3 eV. It is also clear from 
the hgure that the values of f j for CdV204 are always found to be less than ZnV204 and 
MgV204 compounds for whole range of U, which indicates that it is least frustrated among 
these spinels. The reason for the different values of f j for these compounds for whole range 
of U is discussed here as. It is also clear from the Eqn. ( 5 ) that the values of fj depends 
on both the contribution of L to the J as well as the values of f s- The reason for the 
different values of L for these compounds as a function of U has been already discussed 
above. Here, we discuss only for f s- In the tetragonal phase, distorted VOe octahedron 
and V4 tetrahedron have different V-O-V angles and V-V bonds, respectively, which are 
different for these compounds. Because of the different values of these bonds and angles 
from compound to compound (means CdV204 to MgV204 to ZnV204), the values of f s 
is expected to be different for whole range of U studied here. For, example the strength 
of the magnetic interactions (or \Ocw\) are small for CdV204 as compared to MgV204 
and ZnV204 because of the large V-V bonds and V-O-V angles in CdV204 as compared to 
other two compounds. Also, Canosa et al. have proposed that the Tat of these compounds 
also increases with decreasing V-V distances.^ Hence, geometrical frustration is expected to 
different for these compounds because of the different values of |6)(7w|) T^ and L. Here, it is 
interesting to compare the order of calculated values of / j with the calculated values of AM2 
as both are attributed as the degree of geometrical frustration in these compounds. The 
values of fj calculated for U= 3 . 3 - 5.3 eV provide the correct ordering of frustration indices, 
which is consistent with that arises from the AM2. Hence, in this range of U, a long issue 
of the inconsistency about the degree of geometrical frustration has resolved, which indicate 
the importance of orbital degrees of freedom in these compounds. This range of U also 
appear to make sense in the light of the work of Canosa et al., where d electrons of CdV204 
is put in the localized-electron regime; and ZnV204 and MgV204 are put in intermediate 
localization ragime.^ At that point, it is important to note that the calculated values of 
band gap in this range of U are found to be closer to that of experimentally observed values 
for all three compounds. 

Now, in order to evaluate the nearest neighbour exchange coupling constant Jnn between 
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the V atoms of above mentioned compounds, we have considered the FM and AFM (not 
the experimentally observed structure) spin conhgurations. This is because of the fact that 
the experimentally observed sign of exchange coupling constant Jnn gives the information 
about the magnetic ground state of the compound and does not explain it’s spin structure. 
Lowest energy spin structure (i.e. magnetic ground state) is obtained by comparing various 
spin conhgurations of the compound. Hence, the FM and AFM primitive unit cell contains 
four V atoms (V 1 V 2 V 3 V 4 ) with spin tttt Htt (used in the present calculations), 
respectively are sufficient to estimate the Jnn for all three compounds. The atomic and spin 
arrangements of all four V atoms in the AFM primitive unit cell of these spinels are shown 
in the Fig. 5 . Here, we have calculated the total energies of FM and AFM ordering in the 
classical Heisenberg Hamiltonian of the following form® 

H = -JnnY,Si.S, (6) 

i>j 

where. Si and Sj are the spin operators of V sites i and j, respectively with S = 1 for V 
SdJ electrons. In these calculations, we have considered only nearest neighbour interactions. 
In spinels, the magnetic unit cell contains six nearest neighbour interactions among V atoms. 
Hence, the energy of the FM and AFM ordering per unit cell is now expressed as: 


EfM Eq QJnn 


( 7 ) 


and 


EaFM — Eg + 2 J„ 

After solving these two Eqns., we get 


( 8 ) 


Jn.n. 


Eafm — E^ 


FM 
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( 9 ) 


Using Eqn. ( 9 ), we have calculated Jnn for these compounds by varying U from 2 to 
6 eV. The plot of Jnn versus U for these spinels is shown in the Fig. 6. It is evident 
from the hgure that the negative values of exchange coupling constant for these compounds 
indicate the antiferromagnetic interaction for whole range of U. For ZnV204 and MgV204 
the magnitude of Jnn (denoted by \Jnn\) decreases continuously for whole range of U as per 
expectation if The values of \ Jnn\ decreases from ~ 82.4 (~ 90 . 6 ) to ~ 48.7 (~ 30 . 7 ) 
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meV for MgV204 (ZnV204) as U changes from 2-6 eV. For CdV204, \Jnn\ decreases from 
~ 33.4 to ~ 5.6 meV as U changes from 2-6 eV. It is also clear from the hgure that among these 
vanadates, the strength of the exchange interaction is largest for the MgV204 and smallest for 
the CdV204 for U > 2.5 eV. Now, we have estimated the Curie-Weiss temperature {Ocw)s 
for these vanadates from the calculated values of Jnn by using following formula,^ 

{^cw)s = - ^-zJnn ( 10 ) 

3 k B 

where, S'=l (for V 3 (P electrons) and z=6 are the total spin angular momentum and the 
nearest neighbours of V atom (among which the exchange interaction is effective), respec¬ 
tively. Calculated values of {Ocw)s as a function of U for these compounds are given in the 
Fig. 7 . It is clear from the hgure that the absolute values of {Ocw)s [denoted by \{Ocw)s\] 
decreases continuously with increases U from 2 to 6 eV for all compounds. For U < 2.5 eV, 
it is largest for ZnV204 and smallest for CdV204 compound. As U increases from 2.5 to 6 
eV, the order of \{Ocw)s\ becomes MgV204>ZnV204>CdV204, which indicates that the 
strength of the magnetic interaction is highest for MgV204 and lowest for CdV204. Here, 
it is interesting to compare the \{Ocw)s\ with the experimentally observed {Ocw)e^^ for 
3 . 3 < U < 5.3 eV. This is because the consistency about the degree of geometrical frustration 
(which arises from /j and AM2) as well as the the band gap are achieved in this range 
of U as discussed above. Experimentally reported absolute value of {Ocw)exp for ZnV204, 
MgV204 and CdV204 are ~850 K ~600 K and ~400 K, respectively as shown in the Table 
l.SSllt is evident from the hgure that the \{Ocw)s\ is about 2.5 ( 4 . 5 ) times larger than the 
experimental one for ZnV204 (MgV204) compounds for 3 . 3 < U < 5.3 eV. For CdV204, it is 
almost close to the {Ocw)e^-p for this range of U. Here, it is also interesting to see the effect 
of orbital degrees of freedom on the {Ocw)s because the contribution from the orbital part 
of MM is signihcant and can not be neglected for these spinels. Hence, we have replaced the 
total spin angular momentum S by the total angular momentum J in Eqn. ( 10 ), we get 

{^cw)j = ^ 07 ,—( 11 ) 

Skb 

Dividing Eqn. ( 11 ) by Eqn. ( 10 ), we get 

{€>cw)j = ^ 2 — '^~i^cw)s (12) 

for 5 = 1 . 
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Now, Eqn. ( 12 ) is used to calculate the { 0 cw)j) where we have included both orbital 
and spin angular momenta by varying U from 2 to 6 eV. In Fig. 8, we have plotted the 
{Ocw)j for these compounds as a function of U. From hgure, it is clear that the order of 
absolute values of { 0 cw)j [denoted by \{Ocw)j\] is similar to that of \{Ocw)s\ in these 
compounds for whole range of U. For these compounds, the values of \{Ocw)j\ decreases 
continuously as U changes from 2-4 eV. For > 4 eV, the increase in the \ {Ocw)j\ is mainly 
due to the decrease in the absolute values of orbital angular momenta (as discussed above) 
for these spinels. The values of {Ocw)j come closer and closer to the {Ocw)exp for both 
ZnV204 and MgV204 compounds as U varies from 3.3 to 5.3 eV. For ZnV204 (MgV204), 
\{.^cw)j\ increases from ~310 K (~1145 K) to ~560 K (~1550 K) as U changes from 4 to 
5 eV. However for CdV204, it become less than the experimentally reported value, where 
it changes from ~40 K to ~160 K as 1 / increases from 4 to 5 eV. It is well known that 
these spinels belong to the family of geometrically frustrated systems. Hence, the magnetic 
transition temperature (where the magnetic ordering occur) will be different from the Curie- 
Weiss {Ocw)s or {.^cw)j because of the presence degree of geometrical frustration in these 
spinels. 

Now, in order to calculate the magnetic transition temperature for these vanadates, we 
have included the frustration by dividing the \{Ocw)s\ by fs as: 

(T^)s = ( 13 ) 

Js 

for ^=1. 

The plot of {Ti\i)s is shown in the Fig. 9 (a) as a function of U. It is clear from the figure 
that the values of (T]sf)s decreases continuously for all compounds as U increases from 2 
to 6 eV. For ZnV204 (MgV204) compounds, it decreases from ~110 K (~195 K) to ~96 
K (~180 K) when U changes from 3 . 5-5 eV. Similarly, it decreases from ~35 K to ~31 K 
for CdV204 as U varies from 3.5 to 5 eV. From above discussion, it is also clear that the 
values of the (Tj^)s for MgV204>ZnV204>CdV204 as per the experimental results. The 
values of (Tjv)exp for these compounds are also given in the Table 1 . In this range of U, the 
values of (Tn)s are still about 2.5 ( 4 . 5 ) times larger than (TAr)exp for ZnV204 (MgV204) 
compound. However for CdV204, it is close to the experimental one. Here, it is interesting 
to see whether the magnetic transition temperature of these compounds is improved by 
including the orbital degree of freedom in { 0 cw)s- Now, {Ocw)s has become {Ocw)j- The 
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values of magnetic transition temperature (Tn)j are now calculated by dividing \{ 0 cw)j 
by fj as: 


{Tn)j = 


\{Ocw). 

fj 


( 14 ) 


In Fig. 9 (b), we have given the plot of (Tjv)j versus U. From hgure, it is clear that 
the values of (Ta?)/ are also largest for MgV204 and smallest for CdV204 compound for 
whole range of U studied here. Here, it is important to note that the (Tn)j for these 
compounds depend on the ratio of \{Ocw)j\ and fj and hence depends on as well as 
the contribution of L to J. The values of and L depends on the V-V, V -0 bonds 
and V-O-V angles of distorted VOe octahedra, which are different for these compounds as 
discussed above in more detail. Among these spinels, lowest value of (T]y)j for CdV204 
is expected because of large V-V, V -0 bonds and V-O-V angles as compared to MgV204 
and ZnV204 compounds. Also, the order of iT]s[)j remains same as that of {Tn)s and 
(rAr)exp- The values of {Tn)j are ~150 K, ~60 K and ~22 K for MgV204, ZnV204 and 
CdV204, respectively for 3 . 3 < U < 5.3 eV. Now, it is interesting to compare the (Tjv)s, 
{T]sf)j and (TAr)exp for these spinels. It is evident from the Table 1 that the values of 
{T Ar)exp for MgV204, ZnV204 and CdV204 are ~42 K, ~40 K and ~35 K, respectively. For 
3 . 3 < U < 5.3 eV, the (T7v)j is more closer to (T7v)exp as compared to the (Ttv)^ indicates 
the importance of orbital degree of freedom for these compounds. From above discussion, 
it is also clear that the iTisf)j is found to be close to the (T7v)exp for ZnV204 compound. 
However, for CdV204, (Tjv) j is less than the (T N)exp- Here, it is important to note that the 
calculated values of frustration indices for CdV204 compound are very small as U changes 
from 3.3 to 5.3 eV, which indicates that the system is unfrustrated. Hence, it is not necessary 
to include fj for calculating (Tn)j in this system. Now, { 0 cw)j ~ (^iv)j ^ind is found 
to be close to the experimental result. However, only for presenting the order of {Tn)j for 
these compounds, we have included the fj in CdV204. For MgV204, (Tn)j is still about 
3.5 times larger than (Tjv)exp- The various reason for (Tjv) j>(Tjv)exp for this compound 
are discussed below: (i) the magnetic transition temperature as well as the Curie-Weiss 
temperature depends on the purity of the sample. For example, (TAr)exp reported by the 
Nishiguchi et al. is ~45 K for MgV204, which is different from the Mamiya et al. l ^^psps ES] 
(ii) the values of /j (which are used for calculating the (Tn)j) are calculated from the 
Eqns. ( 1 ) and ( 5 ), which shows the dependence of f j on the experimentally reported values 
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of {Ocw)e^p and (T7v)exp- Also, the estimation of {Ocw)exp depends on the linear htting of 
the reciprocal magnetic snsceptibility versns temperatnre data, which is also not consistent 
for this componnd. For example, {Ocw)e:^p estimated by Mamiya et al. is ~600 K, which is 
fonnd to less than as reported by Blasse et al. and (hi) the local fnnctional (LDA+t/) is 
not fonnd to treat the magnetic interactions properly even in MnO, where it overestimates 
the magnitnde of these interactions as compared to experimental resnlts.^ Hence, it is 
expected that the {Ocw)j is greater than the experimental one for above mentioned spinels, 
which may be also responsible for larger valnes of (Tjv)j. At last, we conclnde that even 
having the above mentioned ambignities, the magnetic properties of all three componnds are 
well explained for 3 . 3 < U < 5.3 eV. Hence, the present stndy clearly shows the importance of 
orbital degrees of freedom in nnderstanding the geometrically frnstrated magnetic behavionr 
of these componnds. 


V. CONCLUSIONS 

In conclusion, we have made an attempt to resolve the issue related to the degree of 
geometrical frustration present in AV2O4 (A = Zn, Cd and Mg) compounds that arises from 
the two experimental data (i.e. frustration indices and magnetic moments) by using detailed 
LDA+t/+SOC calculations. The spin-orbit coupling was found to play an important role in 
deciding the magnetic properties of these compounds. The magnitude of maximum calcu¬ 
lated values of orbital magnetic moment per formula unit for CdV204 (~ 1.74 /iB)>ZnV204 
(~ 1.54 /iB)>MgV204 (~ 0.92 ^b)- The difference between the calculated and experimentally 
observed magnetic moment (AM2=Mtotai-Mexp) was observed to be largest for MgV204 and 
smallest for CdV204 compound. The values of the frustration indices (/ j) by including the 
spin and orbital part of magnetic moment were also found to be largest for MgV204 and 
smallest for CdV204 compounds for U= 3 . 3 - 5.3 eV. The consistency about the degree of ge¬ 
ometrical frustration (which arises from AM2 and fj) was obtained for this range of U. The 
order of the absolute values of nearest neighbour exchange coupling constant (Jnn) between 
V spins were found to be MgV204>ZnV204>CdV204, which show that the magnitude of 
Curie-Weiss temperature [{Ocw)j] were highest for MgV204 and smallest for CdV204 for 
3 . 3 < U < 5.3 eV. In this range of t/, the order of magnetic transition temperature [(TAr)j] 
as well as the band gap were found to be in accordance with the experimental results for 
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these compounds. At last, we conclude that all the magnetic properties studied here were 
found to be well explained for 3 . 3 < U < 5.3 eV. 
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TABLE I. Experimentally observed Curie-Weiss temperature [i0cw)exp] (K), magnetic transition 
temperature [(Tjv)ea:p] (K), frustration index (fs), magnetic moment per formula unit [Mexp] (m_b) 


and calculated values of orbital magnetic moment per formula unit [Morbitai] (/^b) for ACr 204 

(A=Zn, Mg and Cd) componnds. 

Compound 

{0Cw)ex^ 


yjlS 

M El 

J-v-i-exp 


ZnV204 

-850 

40 

21.3 

1.26 

-1.50 

MgV204 

-600 

42 

14.3 

0.94 

-0.40 

CdV204 

-400 

35 

11.4 

2.38 

-0.40 
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FIG. 1 . (a) Atomic arrangements of V and O atoms of AV2O4 (A = Zn, Cd and Mg) compounds 
in tetragonal unit cell. Each V atom is located at the center of edge sharing VOg octahedra (V -0 
bond is not connected for the sake of clarity) and (b) a schematic energy level diagram showing 
the splitting of d levels due to crystal field and spin-orbit coupling. 
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FIG. 2. Spin (Mgpin), orbital (Morbitai) and total (Mtotai) magnetic moments per formula unit as 


a function of U obtained from LDA+U+SOC calculations for (a) ZnV204, (b) MgV204 and (c) 


CdV204 compounds. 
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FIG. 3 . AM2 versus U plot for AV2O4 (A = Zn, Cd and Mg) compounds, where AM2=Mtotai-Mexp- 



FIG. 4. The variation of the frustration indices {f j) with respect to U for AV 2 O 4 (A = Zn, Cd 
and Mg) compounds. 
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FIG. 5. The atomic and spin arrangements of all four V atoms in the antiferromagnetic primitive 
unit cell (only used for calculating the nearest neighbour exchange-interaction parameter (Jnn)) of 
AV 2 O 4 (A = Zn, Cd and Mg) compounds. 
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FIG. 6. The nearest neighbour exchange-interaction parameter (Jnn) obtained by calculating the 
total energies per unit cell of ferromagnetic and antiferromagnetic conhgurations onto a classical 
Heisenberg model as a function of U for AV2O4 (A = Zn, Cd and Mg) compounds. 



FIG. 7 . Curie-Weiss temperature [(Ocw)^] estimated from the calculated values of the nearest 
neighbour exchange-interaction parameter (Jnn) versus U for AV2O4 (A = Zn, Cd and Mg) com¬ 
pounds. 
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FIG. 8. The variation of the Curie-Weiss temperature [(0cw)j] as a function of U for AV 2 O 4 (A 
= Zn, Cd and Mg) compounds. 



FIG. 9. (a) The magnetic transition temperatures [(TAr)s'] and (b) [(TAr)j] as a function of U for 
AV 2 O 4 (A = Zn, Cd and Mg) compounds. 
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